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Abstract 

The upper and lower bounds of entanglement of formation are given 
for two mode squeezed thermal state. The bounds are compared with 
other entanglement measure or bounds. The entanglement distillation 
and the relative entropy of entanglement of infinitive squeezed state are 
obtained at the postulation of hashing inequality. 
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Quantum entanglement plays an essential role in all branches of the emerg- 
ing field of quantum information^]. Recently a great deal of attention has been 
devoted to the quantum information processing with canonical continuous vari- 
ables and various protocols for quantum communication and computation are 
developed based on continuous variables |Sj 0]- With continuous variable op- 
tical system these protocols can be implemented. The efficiency of entanglement 
manipulation protocols critically depends on the quality of the entanglement 
that one can generate. For practical purpose [5] one needs to determine the 
entanglement degradation in transmission of two mode squeezed vacuum state 
light through absorbing fibers. It is therefore essential to be able to quantify 
the amount of entanglement in systems with continuous variables. 

Of all the continuous variable quantum states, quantum Gaussian states are 
of great practical importance, since these comprise essentially all the experi- 
mentally realizable continuous variable states, theoretically they might also be 
simple enough to be analytically tractable for the problems involved. Gaussian 
state is completely specified by its mean r\ and its correlation matrix a, where 
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the mean can be dropped by local unitary operation so that is irrelevant for 
entanglement problems. The density operator of Gaussian state can be charac- 
terized [HI by its quantum characteristic function \ ( z ) — exp[i?7 T z — ^z T az] , 
and any Gaussian state of two modes can be transformed into what we called 
the standard form, using local unitary operation only [7] [S] . The correspond- 
ing characteristic function has mean rj = and the correlation matrix has the 
simple form 

61, 0, ci, 
0, bi, 0, c 2 
ci, 0, 6 2 , 
0, c 2 , 0, 6 2 

A particular simple subset of two mode Gaussian states is the two mode 
squeezed thermal states (TMST), with b\ = b 2 = b, c\ = — c 2 = c. The density 
operator of TMST is psr = S% (r) P0S2 (r) , where £2 (r) = exp[— r(af aj" — 
aia^)] is the two mode squeezed operator with squeezing parameter r > 0, and 
po is the product of two thermal states A and B whose average photon numbers 
are Na — Nb = N. The relationship between the two expressions of TMST 
is b = (N + 1/2) cosh 2r, c = (N + l/2)sinh2r. One can use the unit scale 
parameters A = tanhr and v = N/(N + 1), then the inseparability criterion has 
a very simple form of A > i>[7] [Sj. 

Even for the two parameter TMST, however, until now the amount of the 
entanglement in several main measures remains an open question. In this paper, 
I will study the upper and lower bounds for the entanglement of formation of 
TMST. First I will establish the lower bound of entanglement of formation, 
then the upper bound, then compare them with other entanglement bounds 
and measure, at last draw the conclusions. 

Now let us consider the lower bound of entanglement of formation. It was 
proved that entanglement of formation is nonincreasing under local operation 
and classical communication This presents us a nature way to obtain the 
lower bound of the entanglement of formation of a mixed state. By performing 
a special kind of local operation and (or) classical communication to the state, 
converting it to a state whose entanglement of formation is well known, the 
entanglement of formation of the resultant state will provide a lower bound to 
that of the original state. 

For TMST, the simplest way of conversion is to map it on a two-qubit system 
by means of local operation. One of the operation for such purpose is with the 
Hermitian spin one-half operators 10 Sj , p. — A, B for A, B parts respectively, 
SI + iSg = 2 £™ =0 |2m) MM (2m + 1| , ^ = £^ =0 (-1)" 1 \m) ^ (m\ ,which sat- 
isfy the Pauli matrix algebra: [Sf , Sj] — 2ieijk5 l _ l „S£. 

Let us denote the mapped state with p' and assign the following qubit density 
matrix p' A to the original reduced density matrix pa, p'a — \ {Ja + S A ■ a) , where 
S A -a = Y^i=i Tr (paS a ) o~i, Ui are Pauli matrices and Ia is the identity oper- 
ator on the HUbert space of qubit A. 



(1) 
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The map of bipartite density matrix is|l(J| 

p' = ~ ^Ia® Ib + S A ■ a- <g> I B + Ia®S B ■ t ^ (Jl ® CT J j 

, with t l3 = Tr (pSfS?) . 

For TMST density matrix psTi its reduced state p^ = Tr^psT turns out 
to be a one mode thermal state of the form (1 — v') J2m=o v ' m \ m ) ( m l W1 th 
v' = jt£p ■ It is clear that Tr (paSi 2 ) = 0, and a direct calculation shows that 

Tr (paS^) = jt^t- Meanwhile the two mode squeezing operator 62 (r) will gen- 
erate or annihilate the same number of bosons in mode A and mode B, it will not 
change the difference of the occupation numbers of the two modes, one readily 
gets tij = 0, for i ^ j. As for the diagonal elements ta, from the definition of , 
it follows that tn = -*22 = 2Tr £m,„=o Pst |2m, 2n) (2m + 1, 2n + 1| .Making 
use of the coherent state representation of pgr , after some algebra, one obtains 

*33=(^f) 2 and 



00 2 min(m,n) 

/u=2C £ £ 

m,n— Z— 



2n+ 1 
2m + 1 



2to+ 1 
i + 1 



2/? 
J 



2i+l 2m+2n- 



2/ 



with Co = — i_,> A ii 



The qubits map of TMST is 



UJ — — ; TTT" • 

1 — V A 



(2) 



(l+^^l+A^) 1 

0, 

0, 



0, 

V 

0, 

0, 



0. 

0, 

V 

(l+u)*' 

0, 



v 2 +X 2 



(i+vy(i+x 2 ) 



(3) 



The concurrence of p' will be C = max{0,tn — Ivj (1 + v) }. Denote 
x — I (l + Vl — C 2 ) , the entanglement of formation of p 1 reads 



E (p 1 ) — h (x) — — xlogx — (1 — x) log(l — x). 



(4) 



Throughout this paper the base of log is 2. One then obtains the lower bound 
of entanglement of formation : Eif (psr) = E (p r ). 

We now turn to the upper bound for entanglement of formation of quantum 
Gaussian state and specifically TMST. The entanglement of formation is defined 
as follows UJ. Given a density matrix p of bipartite quantum system, consider 
all possible pure state decompositions of p, that is, all ensembles of states |V>i) 
with probabilities pi, such that p = ^2 t Pi \ipi) (ipi\. For each pure state \ip), 
the entanglement E (?/>) is defined as the entropy of the reduced state of \ip}- 
The entanglement of formation of the mixed state p is then defined as the 
average entanglement of the pure states of the decomposition, minimized over 
all decompositions of p : E (p) — xixm^.piE (ipi). 
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From the definition, it is known that any ensemble as a decomposition of p 
has an average entanglement no less than the entanglement of formation of p, 
and so that is a candidate in providing an upper bound for the entanglement of 
formation of p . The aim of us is to construct a series of decompositions, and 
within the series find out the decomposition of minimum average entanglement. 
Because the series is just a subset of the set of all possible decompositions, the 
minimum within the series is a local minimum. It provides a relatively tight 
upper bound for global minimum, the entanglement of formation of p. 

Among all the decompositions, Gaussian state decomposition is mathemat- 
ically easy to deal with. Rather than decompose a Gaussian state to its pure 
Gaussian state ensemble, let us generate mixed Gaussian state from a pure 
Gaussian seed. In the viewpoint of quantum channel, the way of generating 
mixed Gaussian state from a pure Gaussian state is just the way of transmit- 
ting the pure Gaussian state through quantum Gaussian channel ^2- So that 
the problem of how to generate mixed Gaussian state is physical, at least in 
principle. Suppose a pure Gaussian state is given with correlation matrix a s 
and zero mean. Applying local unitary displacement operation on it yields 
another pure Gaussian state, with the same correlation matrix but a nonzero 
mean value or displacement 77. The two states have the same amount of en- 
tanglement because local unitary operation does not change entanglement of 
pure state share by two parts A and B. If the displacement 77 has a Gaussian 
probability distribution, then the mixture of all these Gaussian states will be a 
Gaussian state[o]. Let P be a Gaussian probability distribution with zero mean 
(for convenience) and covariance matrix M for variables (r) A , rf A , rig, rf^) =rj T . 
The quantum characteristic function of the mixed state pp is 

/I 1 
exp(i7y T z — —z T a s z)P (drj) = exp[— —z T (M + a s ) z]. (5) 

The correlation matrix of pp is ctp = M + a 8 . Normalizable of probability 
measure requires that M should be positive semi-definite, M = is included 
to correspond to <5 function type of probability distribution. So that arbitrary 
Gaussian mixed state pp can be decomposed into an ensemble of pure Gaussian 
states as far as the difference of the correlation matrices is positive semi-definite. 
That is 

M = a P - a s > 0. (6) 

So if the separable pure Gaussian state correlation matrix a s is found under the 
restriction of Eq. JfjJ), state pp should be separable. If it is impossible to find 
such an a s , then one has to turn to find an entangled pure Gaussian state cor- 
relation matrix a s to fulfill Eq.©. The minimum of the entanglement of those 
pure entangled Gaussian states will give out the upper bound of entanglement 
of formation of state pp. 

Consider the entangled TMST psr with A > v. Let 
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cosh 2a; , 0, sinh2cc, 

0, cosh2x, 0, — sinh2cc 

sinh2x, 0, cosh 2x, ' ^ ' 

— sinh2a;, 0, cosh 2a; 

be exactly the squeezed vacuum state with squeezing parameter x. M > 
requires that b — | cosh 2a; > 0, (b — i cosh 2a;) 2 — (c — i sinh2a;) 2 > O.This can 
be written as cosh 2 (r — x) < ^_^ 2 .Denote w = tanhw, one gets r — x < w, or 
a; > r — w. The minimum is xo = r — w. Let y = sinh 2 xq, in explicit form of 

( A— u) 2 

A,u, one has y = nZA^TTi^t?! ' r ^^ ie u PP er bound of entanglement of formation 
then will be 




E U f (pst) = 9 (y) = (y + 1) log (y + 1) - y log y (8) 

Let us compare the entanglement bounds obtained above with other bounds 
and entanglement measure. First of them is the upper bound of relative entropy 
of entanglement. The entanglement measure is defined by the minimization 
of the distance of the entangled Gaussian state to the set of separable states 
measured by the relative entropy, its bound was obtained by restricting the 
minimization set to that of separable Gaussian states^- We can further restrict 
it to the set of separable TMST psr — S% (?) p (v) S^" (?) , with A = tanh r < v. 
Then the minimization can be carried out. The upper bound will be 

E ur (Pst) = -2g (N) - 2 log (1 - ?*) - cosh (r - w*) - 1] logv*. (9) 

1 — V 

The minimum reaches at A = v — v* — tanh w* for some v* . 
The second one we should compare with is entanglement measure of loga- 
rithmic negativity for TMST, it reads 

E ln ( P st) = log[(l + A)(l - v)} - log[(l - A)(l + «)]. (10) 

The last one used in comparing is coherent information. The coherent infor- 
mation of bipartite state a with reductions a a and <jb is defined as (a) = 
S (c M ) — S (a) ,for S (cr u ) — S (a) > and I M (a) = otherwise, where S denotes 
the von Neumann entropy. And if for any bipartite state the one-way distill- 
able entanglement is no less than coherent information which is called hashing 
inequality, then one obtains Shannon- like formulas for the capacities JH]- The 
hypothetical hashing inequality reads 

(a) > I B (a) (11) 

where (a) is forward classical communication aided distillable entanglement 
of the state. I do not mean to prove this inequality, but rather to check if it is 
violated for some cases in continuous variable system. It is known that two way 
distillation of entanglement is no less than that of one way|HI, and entanglement 
of formation is no less than relative entropy of entanglement [T?>] . the later is no 
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less than distillation of entanglement 17 . In short one has Ef (a) > E r (a) > 
D<^> (er) > (er). Combining with hashing inequality one then obtains 

E u f (a) > E f (a) > I B (a) , E ur (a) > E r (a) > I B (a) . (12) 

For TMST p S r, the entropy of its reduced state is g (N'), with N' = = 
T, — '^tt-. — rrr. The coherent information is 

(1— u)(l-A J ) 

I B (pst) — max{g (N 1 ) — 2g (N) , 0} (13) 

By comparison of all bounds and measure of TMST I find that: (1) the upper 
bound of entanglement of formation and of relative entropy of entanglement are 
no less than coherent information, so it can be anticipated that hashing inequal- 
ity is also true for TMST or more generally for bipartite Gaussian state; (2) As 
the squeezing parameter tends towards zero, the lower and upper bounds of 
entanglement of formation tend to coincide with each other. (3) As the squeez- 
ing parameter tends towards infinitive, the upper bound of relative entropy of 
entanglement tends to coincide with the coherent information, so if the hashing 
inequality is right, the relative entropy of entanglement, the distillation entan- 
glement and coherent information are all equal for infinitive squeezing state. (4) 
The entanglement measure of logarithmic negativity is at the topmost of all the 
bounds. 

The main points of this paper are as follow: (1) The lower bound of en- 
tanglement of formation of TMST is drawn from the fact that local operation 
can only decreases if not preserves the entanglement, at lower squeezing side, it 
almost coincides with the upper bound. This result leads to the possible deter- 
mination of the entanglement of formation at lower squeezing level. (2)As for 
the upper bound of entanglement of formation of quantum Gaussian state, an 
inequality is derived, for the case of TMST, an analytical expression is given 
out. (3)By carrying out the minimization of the upper bound of relative en- 
tropy of entanglement [S] for TMST I determine the entanglement distillation 
and the relative entropy of entanglement of TMST with infinitive squeezing at 
the postulation of hashing inequality. 
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Figure 1: Comparison of bounds and measure at A — ► 1 side. A = 0.99, dashed 
line for E u f, dashdot for E ur , solid for I B , dotted for Eln 
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Figure 2: Comparison of bounds and measure at A = 0.5, dashed line for E u f, 
dashdot for E ur , the lower solid for Eif, dotted for I B . the topmost solid for 
Eln- 
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Figure 3: Comparison of bounds at A — ► side. A = 0.01, dashed line for E u f, 
dashdot for E ur , solid for Eij, dotted for/ 5 . Entanglement measure Eln is 
omitted for significantly greater than the bounds. 



8 



